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ABSTRACT 

We present a new resummation formula for the Drell-Yan cross section. The formal 
resummation of threshold corrections in Drell-Yan hard-scattering functions produces an 
exponent with singularities from the infrared pole of the QCD running coupling. Our 
reformulation treats such 'infrared renormalons' by a principal value prescription, analo- 
gous to a modified Borel transform. The resulting expression includes all large threshold 
corrections to the hard scattering function as an asymptotic series in a s , but is a finite 
function of Q 2 . We find that the ambiguities of the resummed perturbation theory imply 
the presence of higher twist corrections to quark-antiquark hard-scattering functions that 
begin at ^-qcd/Q- This suggests an important role for higher twist in the phenomonolgy 
of hadron-hadron inclusive cross sections. We also discuss the numerical evaluation of the 
exponent and its asymptotic perturbation series for representative values of Q 2 . 
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1. Introduction 



It has been known for some time that perturbative QCD corrections to the inclu- 
sive Drell-Yan and other hadron-hadron hard scattering cross sections are numerically 
important, even at moderately high energies' 11 . Considerable effort has been devoted to 
developing all-order resummation techniques in order both to study the convergence of 
the QCD perturbation series and to control the size of perturbative corrections to these 

[21 

processes. In a recent paper , we presented an analysis of the resummed threshold correc- 
tions in the Drell-Yan cross section, deriving an explicit expression for the hard scattering 
function directly in momentum space. This new resummation formula organizes all large 
and order unity threshold corrections in the entire region where perturbation theory gives 
the dominant contribution. Because the resummed expression involves integrals over the 
scales of running couplings, it is undefined when these scales reach Kqcd- The purpose of 
the present work is to address these and related issues. 

Beyond specific applications to dilepton cross sections, we are concerned here with 
some very general problems in the application of perturbative QCD to hard processes. A 
complete picture of such processes requires the inclusion of power-suppressed, or 'higher 
twist' contributions, but to include such contributions, we must first define the 'leading- 
twist', or perturbative series' 31 ' 41 . To do so fully is a formidable task, beyond our present 
abilities. The class of large perturbative threshold corrections to dilepton (and other) in- 
clusive hard scattering cross sections, however, is known to all orders. It is, as expected, 
ambiguous. Our goal, then, will be to 'make room' for higher twist by defining the re- 
summation of large corrections at leading twist. At the same time, we may also ask why, 
despite all these large corrections and the ambiguities in perturbation theory, perturbative 
corrections at the one loop level are not totally wrong. While we cannot fully answer this 
question either, we will get some hints toward its resolution. 

The situation for dilepton production hard-scattering functions is analogous, but not 
identical, to the situation for the e + e~ annihilation total cross section 131 141 . In e + e~, the 
ambiguities of perturbation theory may be identified with the gluon condensate (0|F 2 |0), 
which appears in the operator product expansion for this process. Such contributions 
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are suppressed by the power Q~ 4 relative to fixed-order perturbation theory. In contrast, 
we shall identify below the first ambiguity in the resummed perturbative cross section 
at the much larger level of Q" 1 relative to leading power. This suggests an important 
phenomonological role for higher twist in hadron-hadron scattering, at least at moderate 
energies. Another difference between dilepton production and the e + e~ annihilation cross 
section is that in the latter case low order corrections are much smaller. It is possible that 
the occurrence of nonperturbative corrections at lower twist in hadron-hadron scattering 
is related to the larger size of its perturbative corrections. 

Let us now review the basic results of the resummation program applied to the dilepton 
cross section 

hi(pi) + h 2 (p 2 )^U(Q») + X, (1.1) 

with a produced lepton pair of momentum Q^. In perturbation theory, the factorized form 
of this cross section is given by 



da 
W 2 



= a °Yl [ ^ o/hi^o) Q 2 )(Pb/h2 ( x b, Q 2 )uj ab (r/x a x b , a s (Q 2 )) , (1.2) 

J X a Xjj 
ab q 



where ao is the Born cross section, (pa/^ (x a , Q 2 ) is the distribution function of parton a 
in hadron h\, and uj a \, is a short-distance function, or 'hard part'. In eq. (1.2) we denote 
t = Q 2 1 's, with s = (pi +P2) 2 - We emphasize that when we speak of 'leading' and 'higher 
twist' in this paper, we are referring to the hard parts, and not to the cross section as a 
whole, which depends on the interplay of the short-distance functions with the evolving 
distributions. 

In eq. (1.2) the hard parts cu a fr, calculable in perturbation theory, as well as the non- 
perturbative parton densities, are not unique. We shall identify the distribution a // li with 
the contribution, F a i} H , of parton a to a structure function in deeply inelastic scattering of a 
nucleon hi. Then the hard parts uj ab , which directly give a prediction for the normalization 
of the Drell-Yan cross section in terms of the observables F a /f li , contain large corrections. 
These come from large numerical coefficients of 5-functions, and from 'plus distributions'. 
The organization of such large corrections to all orders in perturbation theory, is the 
content of various resummation techniques' 01 161 171 181 . It is convenient, in order to study the 
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large perturbative contributions in eq. (1.2) which come from the region z = r /x a Xb — > 1 of 
the integrals, to take moments and study the large-n region. The only quantities that are 
singular in the z — > 1 limit are the diagonal-flavor quark-antiquark hard parts, for which 
we obtain the exponentiated form 



&qq(n, Q 2 ) = e 2 A(a s (Q 2 ))I(n, Q< 



(1.3) 



where the terms that are singular in the limit n — > oo are contained in 



I(n, Q 2 ) = exp 



1 , X 

, x n ~ x — 1\ f f dy 
— ax 



1 — x 




-Jt- gi (a s [(l-x)(l-y)Q z }) 

y 



+ g2 (a s [(l-x)Q 2 ]) 



(1.4) 



= exp[E(n,g 2 )] 



In eq. (1.3), A(a s (<5 2 )) represents 5(1 — z) contributions, including the exponentiated 
Sudakov 7t 2 terms' 81 . I(n, Q 2 ) contains the exponentiation of all plus-distributions, which 
are the source of growth with n. In ref. 2, we employed an asymptotic expansion for the 
exponent of eq. (1.4), to obtain, directly in momentum space, the resummed hard parts in 
the form 



I(z, a s ) = 5(1 - z) - 



e E (ih> a *) / / l u / / l 



(1.5) 



where 



d 

Pi(n, a s ) = — E(n, a s ) 

omn 



(1.6) 



The central issue that we did not address in ref. 2 was the nature of this asymptotic 
approximation and its relation to non-perturbative effects. This is the main subject of the 
present work. 

Even though u q q is finite at any finite order of perturbation theory, the exponentiated 
form, eqs. (1.4), (1.5), suffers from singularities associated with the behavior of the running 
coupling at small energy scales. Therefore, even though the above resummation formulas 
are useful in reproducing finite-order results' 91 , a direct comparison with phenomenology 
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can be made only if a regularization of soft-gluon effects is supplied. Hence, a sensitivity of 
the resummed predictions to the treatment of this 'soft' region, x — > 1 in eq. (1.4), is to be 
expected 161 . For an infrared safe function, such regions give contributions that are finite 
order-by-order in perturbation theory, but which diverge in the sum to all orders. This 
behavior is, in fact, ubiquitous in QCD 131 [101 . In addition to its applications to the Drell- 
Yan and related cross sections, the resummation (1.5) may thus serve as a 'laboratory' to 
test our ideas about the asymptotic nature of the QCD perturbation series. 

In this paper, we shall propose a finite definition of the exponent E(n, Q 2 ) in eq. (1.4) 
as a principal value integral in the moment variable x. In this fashion we will define 
unambiguously - albeit not uniquely - the perturbative content of large corrections to 
the hard scattering function. We will see as well that any such definition requires that 
nonperturbative corrections to the hard scattering functions begin at order a full 

power of momentum transfer larger than in spin-averaged deeply inelastic scattering, and 
three powers larger than in the total e + e~ annihilation cross section. 

Let us reemphasize that our choice for a principal value resummation is not unique. The 
choice to be given below is closely linked to the specific form for the resummed hard part 
given in eq. (1.4). Other expressions for u) qq , which sum the same sets of large perturbative 
contributions, but which differ in nonleading terms, are possible. For such expressions, a 
different use of the principal value, or even the use of some other prescription, may be 
advantageous. Differences in these prescriptions will be reflected in differences in higher 
twist corrections' 3 " 41 . In fact, there is a whole class of resummation formulas, which can 
be constructed by exploiting the analytic structure of the perturbative running coupling. 
In addition to studying our specific construction, we would like to argue that these issues 
bear further investigation in hadron-hadron scattering cross sections. 

In section 2 we define our principal value resummation formula for oj qq (n : Q 2 ) : and 
discuss its relation to the Borel transform. We relate the exponentiated principal value 
prescription to an analogous prescription for the evolution equation satisfied by u> qq (z, Q 2 ). 
We go on to evaluate analytically the leading (one-loop in the gi and the running coupling) 
exponent in eq. (1.4). We then discuss the consequences of the principal value prescrip- 
tion for higher twist. Here we identify contributions at order 1/Q, associated with the 
ambiguities of the resummed perturbation theory. In section 3, we construct an explicit 
asymptotic perturbation series for the exponent, and in section 4 we explore the exponent 
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and its asymptotic series numerically, again approximating it by the one-loop terms of the 
functions g{. We find that a truncated perturbative series can give an excellent approxi- 
mation to the exponent over a wide range in the variable n, but that such approximations 
require an energy-dependent maximum order, which differs for different terms in E(n, Q 2 ). 
We shall also see that the full principal value exponent increases to a moderate value as n 
increases, and then turns over and decreases as n — > oo. Finally, in section 5 we summa- 
rize our results and conclusions, and discuss the outlook for applying this method to the 
calculation of other hard hadronic processes. 



2. The Principal Value Exponent 

As explained in the introduction, we will construct a resummation formula analogous 
to eq. (1.4), that is calculable without reference to explicit IR cutoffs. We shall deal directly 
with the exponent E(n, Q 2 ) given by (1.4). This defines the resummation in both moment 
and momentum space. 

2.1 Definition 

Denoting the expression in the curly brackets of eq. (1.4) by 

X 

T q<1 (l-x,Q 2 ) = f ^- gi (a s [(l - x)(l - y)Q 2 }) + g 2 (a s [(l - x)Q 2 }) , (2.1) 



we define the principal value exponent for the quark cross section in moment space through 

///■n—l _ i , 
w )r w -(i-C,Q 2 ) 

p 

l+ie ^ 1—ie ^ 

= ~{ J dc(^—±)r q ,(i-CQ 2 )+ J d c(^)r w (i-C,Q 2 )}. 

Q+ie 0-ie 

(2.2) 

Wherever T q q(l — ( 7 Q 2 ) is an analytic function between and 1, we can deform the two 
contours back to the real axis. This is the case at any finite order in perturbation theory. 
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The resummed exponent, on the other hand, contains IR divergences near x = 1 for any 
value of n from the running of the coupling. In this region, the quark cross section defined 
through E(n, Q 2 ) remains well defined with the principal value prescription, although, of 
course, contributions from this region are still infrared-sensitive. 

How may we motivate the prinicpal value prescription? In some sense, the choice is 
quite arbitrary. But some choice is necessary if we are to 'make room' for higher- twist 
effects in physical situations where, as in dileption cross sections, perturbative corrections 
are large. In our opinion, it is precisely such quantities from which we may eventually learn 
the most about the interplay of perturbative and nonperturbative effects in QCD. With 
these observations in mind, we may offer a formulation of principal value resummation 
that is somewhat more general than the specific integrals in eq. (1.4). 

The principal value prescription may be used to define a resummed series whenever the 
summation of a set of perturbative contributions can be expressed as an integral over the 
running coupling. Infrared renormalons' 101 are the simplest examples of this form. This 
set of contributions is identified by a behavior from individual diagrams at nth order of 
a ™ (Q 2 )^2 n '> with Q 2 any fixed scale. Such a series of terms may be generated by expanding 
the running coupling a s (k 2 ) inside the simple integral 

Q 2 

I(a s (Q 2 )) = J alk 2 k 2 a s (k 2 ) (2.3) 
o 

as a power series in a a (Q 2 ) : using the one-loop expression, 

a(k 2 ) = asiQ2) = - (2 4) 

s[ ] l + (6 2 /7r)a s (Q 2 )ln(fc 2 /Q 2 ) &2ln(fc 2 /A 2 )" 1 ' ] 

Many of the properties of infrared renormalons are brought out very clearly by reex- 
pressing the series as a Borel transform, 

ioo 

1(b) = J d(l/a s (Q 2 ))e^ b / a ^I(a s (Q 2 )), (2.5) 
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whose formal inverse is 

00 

I(a s (Q 2 )) = J dbe- b / a ^I(b) . (2.6) 
o 

The integral I(a s (Q 2 )) is undefined as it stands, because of the singularity in the running 
coupling a s (k 2 ) at k 2 = A 2 . This singularity is directly reflected in the n\ behavior of the 
expansion in a s (Q 2 ), and in a consequent singularity on the positive real axis in the Borel 
transform, 1(b), at b = 27r/&2- In Ref. 3, it was shown that, given what we know now 
about the perturbative series, this is the leftmost singularity of the transform. Assuming 
the inverse transform, eq. (2.6), has been defined, such a singularity may be expected to 
contribute only at the level to the cross section, corresponding exactly to the contri- 
bution of the gluon condensate, (0| F^F^ |0). One way to define the perturbative, and 
hence nonperturbative, contribution at this level is to define the inverse Borel transform 
as a principal value' 31 . In fact, for the one-loop running coupling the resulting expression 
is exactly what we find by treating the original integral in eq. (2.3) as a principal value in 
k 2 , and by changing variables to 

b' = ln(k 2 /Q 2 ) . (2.7) 



A similar change allows us to reproduce two-loop results as well. 

In the somewhat more complex situation of eq. (1.4), integrals of the running coupling 
appear in an exponent, and the relation of a principal value in the moment integrals to an 
inverse Borel transform for I(n, Q 2 ) is not quite so simple. The prescription can still be 
given a reasonably general motivation, however, in an evolution equation satisfied by u) qq -. 



Up to corrections that are finite in the z — > 1 limit, the dilepton hard-scattering 



function u> qq -(z, Q 2 ) obeys an evolution equation of the form 



[5] 



d 



01n(l/(l-z)) 



^)^ + l]o/«(z,Q 2 ) 



l 

-fdy 

Z 



(i-y) 



+ 



Uqq(y ~ Z, Q 2 ) + g2(a s (Q 2 ))uJqq(z, Q 2 ) 



(2.8) 
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whose solution is 



qq 



oo n „ ^ r 





^(a s ((l-a:,)(l-^)Q 2 ))' 



n=0 i=0 ' 



1 - Hi 



+ 



(2.9) 



+ 



g 2 (a s ((l - Xi )Q 2 )) 
(!-.'•;) 



fe=i 



Moments of this expression with respect to n yield, up to corrections of order 1/n, 
uj q q{n,Q 2 ) in eqs. (1.3) and (1.4). To verify that eq. (2.9) satisfies (2.8), we may con- 
vert the derivative of the delta function with respect to ln(l — z) into a sum of derivatives 
with respect to the x^, and then integrate by parts. The remainder of the reasoning only 
requires the definitions of the running coupling and the plus distribution. 

In terms of the formulation for u q q just given, we see that the singularity in the running 
coupling is already present in the evolution equation, (2.8). If we choose to define the y 
integral of this equation as a principal value, its solution (2.9) is then given as a sum of 
products of principal value X{ integrals, and its moments become exponentials of principal 
value integrals. (The delta functions may be considered as referring to the real parts 
of the X{.) We may, therefore, wish to consider our principal value prescription for the 
exponentiated resummation as grounded in a principal value definition of the integro- 
differential equation whose solution it is. 

Some basic properties of the principal value resummation formula are immediately 
obvious. First, the exponentiated cross section in moment space, eq. (1.4), is real (for real 
n). This is because the contour P in eq. (2.2) is a sum of two mirror-symmetric contours 
with respect to the real axis, and the integrand is real on some portion of the real axis. 
Also, E in eq. (1.4) is finite by the definition of T qq -, eq. (2.1). 

At this point, we may come back to the Borel transformation . All large perturbative 
corrections' 21 are contained in those pieces of the function T qq {l — £, Q 2 ), eq. (2.1), whose 
dependence on Q and Q 2 is of the form 

2 

r«(i - C, Q 2 ) = £ « fc r^ (a in[i/(i - c)]) , (2.10) 

fc=0 
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where we define 



a = a s (Q 2 )/iv, 



(2.11) 



and where each is an infinite series in its arugment. In our exponent 

£(n,Q 2 ) = -J> fc f rfC(^^)r^(aln[l/(l-C)]) , (2.12) 

fc=o p ^ 

we now expand the power of £, 

c- 1 = (i-a-or 1 = £(i-n)...(i-n+m-i)-(i-cr = E ^-dn^-cr 



m! 14 — ' m! 

m—(j m=0 



(2.13) 

where 'Pochhammer's symbol' is defined as (a) m = T(a + m)/T(a). Performing the 
change of variables 

w = maln[l/(l-()} (2.14) 
we arrive at the following form for the exponent: 

i 2 °° /i _ \ 

E ( n ,Q 2) = _I^ Q > £ (Lz^= 7 w (Q;m) . (2 . 15 ) 

k=0 m=l 

Here, 7^ is given in terms of T^, eq. (2.10), as 

-y$(a; m) = J dwe~ w ' a T^(w/m) , (2.16) 
P' 

(k) 

where the principal value contour P runs between and 00. This definition of 7^ exactly 
coincides with the inverse Borel transform, defined as a principal value 

2.2 Evaluation of the Exponent 

We are now ready to discuss the explicit evaluation of the function E{n, Q 2 ), eq. (2.2), 
which is defined in terms of the functions gi in (2.1). The functions gi are simple expansions 
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in terms of the running coupling constant ct s [AQ 2 ], with A a scale, 



3=1 



The lowest order numerical coefficients, that determine the large perturbative behavior, 

, [5] [6] 



are 



9 <« = 2C F , 4" = ~C F , 9 < 2 » = C F [C A (% - \) - *%] , (2.18) 

where n f is the number of flavors. We can relate the resummed expression to one contain- 
ing the fixed coupling a s = a s (Q 2 ), by reexpressing the running coupling at the scaling 
variables (1 — — y) and (1 — Q in terms of the fixed one, through the beta function. 
As remarked in ref. 2, all the large perturbative corrections in the resummation, of order 
unity or greater in the range a s ln[l/(l — z)\ < 1, are contained in terms that have at least 
as many powers of the logarithm of the momentum scale as of the fixed coupling. This 
means, in turn, that we only need to keep terms that are leading and next-to-leading in 
the gi and the QCD beta-function. We now construct these terms explicitly. 

Consider the renormalization-group equation for the running coupling a = a s /iv, 

-a[XQ 2 ] = -b 2 (a[XQ 2 ]) 2 - b 3 (a[XQ 2 ]f (2.19) 



<91n(A) 



where, in QCD, 



b 2 = (33-2n/)/12 , 6 3 = (306 - 38n/)/48 . (2.20) 
The solution of eq. (2.19), with the initial condition q;(A)|a=i = ot = a s {Q 2 )/-K, is 



-l 

a(X)/a= l+ab 2 \n(X)-(ab 3 /b2)[Ha(X))/a)-\n((b 2 /b 3 +a(X))/(b 2 /b 3 +a)) . (2.21) 



We can solve this transcendental equation iteratively, keeping only leading (a k ln k A) and 



11 



next-to-leading (a k ln k 1 A) powers. These terms are given by 

/ i 1 , , ,, \ ln(l + abo In A) 

«A/a = - — — - (ab 3 b 2 Hr r f ,J • 2 - 22 

v " 1 + 0*2 hi A v ' ' (l + a6 2 lnA) 2 v ' 

Combining eqs. (2.17) and (2.22) we find 

M) . f 9 f aj(T _L_ _ {ah/ ^<L^ r . ( , 23) 

The leading and next-to-leading terms, which give all large perturbative corrections, will 
come from the j ' = 1, 2 terms only. Keeping just the leading and next-to- leading pieces in 
these terms, we can write gi ~ gf + gf L with 

sf ( a[ AQ^])= 9 f) QTT ^ (2.24) 



sfW ]) - ^gig.^^ • (-5) 



Note that the distinction between 'L' and 'iVX' refers only to the loop order in and not 
to the powers of logarithms of n in E(n,Q 2 ). In fact, g^\ for instance, gives the same 
maximum power of logarithms of n as g£\ For now, however, it will be instructive, and 
simpler, to focus on the one-loop terms only. 

Let us now proceed with the analytic evaluation of the exponent 

E(n,a) = - j '_~ ^ ^-(l - C, Q 2 ) ^ E(n, a) L + E(n, a) NL . (2.26) 
P 

Here and below, we replace Q 2 by a as the second argument of E. El, as specified by 
(2.1) and (2.24), is 

E(n,a) L = a(g[ 1) I 1 -g { 2 1) I 2 ), (2.27) 

with 
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and 



//~Ti 1 1 

<%{ \_~ C ) ln ( X + (V*) ln (! " 0) , (2-29) 

P 

6(t)5 / JC (trlzl) i 1 (2.30) 
p 



Here we define 



* = 1/(0*2) = MQ 2 /A 2 ) , (2-31) 
where we have used the one-loop running coupling, as in eq. (2.24). 

The £ integrals in I\ that result from the expansion (2.13) are readily carried out by a 
change of variables to ln(l — £), followed by integration by parts, which gives 

h=y { ™ } ™ {£{mt) - 2S(mt/2)} , (2.32) 
z — ' m!m z 

m=l 

where we define the combination 

£{x) = xe- x Ei(x) , x > , (2.33) 
with the Exponential Integral defined as the principal value integral 1 " 1 

X 

/ y 
dy— . (2.34) 
y 

— oo 

Similarly we find 

h=y (1 ~, n)m g(mt) . (2.35) 

m=l 

Hence the leading exponent can be written as 
E(n, a) L = - ag U £ [2S(mt/2) - £{mt)] - ag? £ ^^e(rnt) . (2.36) 

m=l m=l 

Notice in the above formulas that the moment dependence is entirely contained within the 
Pochhammer symbol, and is defined for complex n. The ct s ((3 2 )-dependence, on the other 
hand, is contained within the functions £(mt), which are defined for any value of t (a s ), 
however small (large). 
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2.3 Higher Twist in the Resummed Exponent 



Eq. (2.36) affords a direct estimate of the ambiguities implicit in the pertrubative series. 
They are, as expected, of higher twist. We have treated these ambiguities by the prinicpal 
value prescription in eq. (2.34). Any other prescription for defining the integrals would 
differ in the treatment of the singularity at y = 0. But the integral at y = is proportional 
to an exponential of —t, and is hence suppressed by a power of A/Q. The relevant powers 
may simply be read off from the arguments of the S in eq. (2.36). The minimum supression 
is from S(t/2) in the m = 1 term, which behaves as (A/Q) ln(Q 2 /A 2 ). At this first 
nonleading power, perturbation theory is already ambiguous, and nonperturbative effects 
must come into play. 

An alternative to the principal value prescription is to simply cut off the £ integral in 
eq. (2.1) at a value large enough to avoid the singularities of the running coupling. It is 
perhaps worthwile to illustrate the relationship between these two approaches. In ref. 6. 
it was shown that the IR cutoff dependence in the hard part, despite being numerically 
significant, is higher-twist, although the precise powers were not determined. These powers 
may be easily determined if we take eq. (2.36) as a starting point, redefining the integrals 
that define the £'s to reflect the cut-off. 

For example, consider an exponent defined by cutting off the ( integral at some min- 
imum value of 1 — ( = £. For the integral I2, for example, the corresponding regulated 
expression would be 

geg) = g (1 ~ w) m j ( reg)m (2 3?) 



m 



with 



mln(l/&>) 



s / «-' TTTI ^ = *® / . (2.38) 

imin mln(£ m4n / '£0) 

In the second form, £0 = exp(— t) is the position of the pole in £. Since £ m i n is designed 
to make the integral finite, it must be chosen to be larger than the location of the pole. 
To separate only the nonperturbative region, however, we take it to be the same order of 
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magnitude as £o- We parameterize the relation by 

a = ln(£ min /£ ) , (2-39) 
with < a < 1. Then, using eqs. (2.35) and (2.38) we obtain: 

ma 

/x 
dx^— = te- mt Ei(ma) . (2.40) 

— oo 

For a given a, the difference between the cutoff integral and the principal value integral is 
exponentially suppressed by the same function of the fixed coupling constant that appears 
in the principal value expression (i.e., by the m-th power of A 2 /Q 2 ). Exactly analogous 
results hold for I\. 



3. The Asymptotic Series 

Even though eq. (2.36) gives the result for the exponent in a form appropriate for 
numerical evaluation, it is also of interest to study the asymptotic expansion of E in 
a s (Q 2 ). Notice that the exponential integrals in eq. (2.36) have a perfectly well-defined 
Taylor expansion [n] 

Ei(mt) = ^ + \n(mt) + ^2^^-, (3.1) 

n=l 

which actually converges better than an exponential for any value of mt. On the other 
hand, such a Taylor expansion would reproduce an infinite series of inverse powers of a s . 
The issue we would like to address here is how to recover a perturbative series for the 
exponent, eq. (2.36). 

We can obtain for the special function £(x), eq. (2.33), after repeatedly integrating by 
parts, the asymptotic expression 

E(mt)~Y-?— , (3.2) 

with N chosen to optimize the approximation. In this asymptotic series, we see explicitly 
the pi behavior at pth order, characteristic of infrared renormalons [101 . Their presence is 
a direct consequence of the singularity in the perturbative running coupling. 
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As we shall show in section 4, a common optimum N for all m can be determined 
numerically for each of the three sums in eq. (2.36). However, the three optimum numbers 
of asymptotic terms differ for the three sums, so we will use the (clumsy but we hope clear) 
notation N = {N[2I(t/2)], N[I(t)], N[I 2 (t)}} for the three sums. Using the asymptotic 
expansion eq. (3.2) in the expression for El, eq. (2.36), the leading exponent becomes a 
finite perturbative sum: 



E(„,a,N) L = -a 9 {»{ £ ^'"WI'Ey 

^ p=0 m=l 
N[I(t)] oo n s v N[I 2 (t)] 

- E E 

p=0 m=l p=0 m=l 

(3.3) 



We can reexpress the infinite series in the moment variable in terms of a plus-distribution, 
through the identity 



T (l-n) m = (-1)" } dxx n-i( ^--) ] (34) 
^ m\mP+ 1 T(p + l)J dXX V l-x )' { ' 

171=1 n + 







These integrals can in turn be expressed, for large n, as polynomials in Inn. Using the 
relation 

d \P, 



we obtain 



ln p (l-x)= lim ( — ) (l-x) e , (3.5) 

e^0+ \0€/ 



£ (1 ,""!? = lim J" 1 )' (|y{B(n,e)-i|, (3.6) 
m l m p+i e ^o+ T(p+ 1) V^e/ I v ' ; eJ' v ; 

m=l 



where B is the Beta function. Using Stirling's formula for large n, we find 

B(n,e) ~e- £lnn r(e) , (3.7) 



and we can approximate the infinite series in the moment variable, n, by a polynomial in 
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Inn: 

t^t^-^^tfwn. (3.8) 

m=l j=0 J 

The numbers Ck above are the standard coefficients of the Taylor expansion [n] , 



oo 

r(l + z) = J2ckZ k - (3-9) 

k=0 



To summarize, the leading exponent may be approximated by a finite perturbative 
sum E(n, a, N) x, 

E(n,a) L = E(n,a,N) L + A(n,a,N) L , (3.10) 

where the second term on the right is a remainder which, since it is not expressible in 
terms of powers of ct s , contains a higher- twist contribution. We shall discuss the numerical 
size of the remainder in section 4. 

The truncated expansion E(n, a, N)l defines a convergent, resummed perturbative 
series for the exponent E, and hence for the hard-scattering function. N is defined in 
eq. (3.10), to minimize Ax,. Increasing iV further will result in a smaller accuracy (larger 
A/,). Inversely, for a given accuracy, there is a maximum a s beyond which the approxima- 
tion of eq. (3.10) breaks down for any N, i.e., the remainder is outside the desired accuracy. 
These features of the exponent represent a 'nonperturbative barrier' in the accuracy of the 
asymptotic approximation, which is what we might expect. 

The leading asymptotic exponent is given by eq. (3.8), in (3.3), as 



JV[2/(t/2)]+l p+l N[I(t)] + l p+l 

E(n,a,N) L = £ a^4[2/(t/2)]Wn + £ a'^^Wn 

N[h{t)+1] p V ' 

+ E « P I>jyj 2 (t)]ln'n 
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with 

sl p [2I{t/2)\ = -g[%^\-iy+i { -Pj^2% +1 _ j , 

sfp[i(t)] = g{%r\-iy +j{ - E jr 1 c p+1 - J , (3.12) 

Again, the A's to be kept in the asymptotic expansion are fixed by minimizing the re- 
mainder in eq. (3.2) for a fixed value of a s (Q 2 ). We may in principle use eqs. (3.11) and 
(3.12) (along with the corresponding results for E^l) to approximate the resummed nor- 
malization by an analytic expression . Alternately, we may use the full principal-value 
exponents. Using either the exact E, or its representation as an asymptotic series, the 
large perturbative corrections to the hard scattering function in momentum space can 
be found from eqs. (1.5), (1.6). 



4. Behavior of the Exponent 

In this section we shall address several important numerical issues. Using the 'leading' 
exponent E(n, a)i, we will illustrate the behavior of the principal value definition and its 
approximation as an asymptotic series, E(n, a, N)i. We defer a similar discussion for the 
next-to-leading exponent, E(n, a, N') _/vx, as well as numerical results for cross sections, to 
future work. We will see that the principal value exponent behaves in a relatively mild 
fashion for all n, over a wide range of Q, and that the asymptotic series can give a good 
approximation to it unless n is very large. 

In the following, we present numerical results for three values of Q, taking A = 
200MeV. To be specific, we have used the values 

a(Q = 5GeV) = 0.075 , a(Q = lOGeV) = 0.061 , a(Q = 90GeV) = 0.039 

[21 

along with 

b 2 (n f = 4) = 2.08333 , = 8/3 , = -2 . 
The optimum A, determined numerically, are shown in table 1. 
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Several details concerning the numerical calculations may be found in the appendix. 
Notice from eqs. (3.11), (3.12) that, for a given power p of the coupling constant, the 
maximum power of In n in each of the three sums has coefficients 



s 



P(P+1) ' 

*i, P m)]=gf h 



,L lT /,m _ „(1)&2 1 



TABLE 1 

Optimum numbers of asymptotic terms as a function of n and a 



a 


n 


N[2I(t/2)} 


iV[J(t)] 


NW)] 


0.075 


20 


1 


5 


5 


55 


40 


1 


•5 


5 




60 


1 


•5 


5 


0.061 


20 


2 


6 


6 


55 


40 


2 


6 


6 


55 


60 


2 


6 


7 


0.039 


20 


4 


13 


11 


55 


40 


4 


12 


11 


55 


60 


4 


9 
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From table 1 we see that the leading power of Inn comes by far from the N[I(t)} and the 
N[l2(t)] terms, and the corresponding coefficients, as shown in eq. (4.1) are both negative. 
This shows that the leading exponent, expressed perturbatively as an asymptotic approx- 
imation, becomes negative for sufficiently large values of the moment n. In momentum 
space this shows that, within perturbation theory, the corresponding cross section tends 
to a finite limit as z — > 1, since the hard part is the exponential of a quantity that diverges 
to minus infinity in that region of phase space. Of course, our asymptotic approximation 
is not valid at the edge of phase space. The preceding discussion, however, serves to point 
out that the corresponding resummed perturbative cross section is finite in that range, and 
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in no need of extra IR cut-offs. What is more important, we can reproduce this satisfac- 
tory feature of the hard part without even restricting ourselves to the perturbative regime. 
Indeed, we may numerically calculate the leading exponent either from eqs. (2.27)-(2.30), 
or from their analytical equivalent, eq. (2.36). The calculation is carried out in moment 
space, but conclusions can be reached in momentum space with the correspondence 

n <-> — . (4.2) 
1 — z 

The calculation of E(n,a)i from the series in eq. (2.36) is pretty sensitive to the 
accuracy with which intermediate quantities are calculated, because there are large can- 
celations involved. Alternately, as we explain in the appendix, it is convenient to evaluate 
the integrals in eqs. (2.28)-(2.30) numerically on a deformed version of the principal value 
contour, which is much more stable from a numerical point of view. This contour, denoted 
by P, is shown in fig. 1. Some details of the integration may be found in the appendix. 
Using these integrals in eq. (2.27), we find the solid curves of fig. 2. On the other hand, 
we may use the asymptotic expression, eqs. (3.11), (3.12), together with values of N like 
those in table 1, to obtain a perturbative asymptotic series for the leading exponent. The 
corresponding values of E(n, a, N)i as a function of n are shown by the dotted curves of 
fig. 2. 

The first thing we note is that the exponent is bounded at fixed Q 2 , and that it reaches 
its maximum at rather large values of n. Beyond the maximum, it decreases monotonically 
toward — oo, due to the dominance in E of the integral h(t), eq. (2.29), which behaves 
as —t ln(2) ln(n) when ln(n) 3> t. Comparison between the solid and dotted curves will 
establish the range of moments n where the perturbative expression is valid, as a function 
of the fixed coupling constant. Below its maximum, E(n, a) l is nicely approximated by its 
asymptotic series, also shown in fig. 2. Beyond this, higher twist takes over and damps the 
exponent, eventually making it negative. The asymptotic series behaves in an analogous 
manner, but decreases only for much larger n. In either case, the limit of the exponent is 
minus infinity as z — > 1, and the integral over z in the cross section is actually finite for 
both the full principal value exponent and its asymptotic expansion constructed as above. 

Introducing the notation 

n\ = n(E(n, a)i = max) , 112 = n(E(n, a, N)i = max J 
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and, with Al given by eq. (3.10), 



fr(„, a)s ' A (?°-?>*l xl00%, 
E(n, a) L 

we may summarize the main features of this comparison in the following table: 

TABLE 2 

The exact exponent versus its asymptotic approximation 



a 


ni 


E(m,a) L 


E(n h a,N) L 


S L (ni,a) 


n 2 


0.075 


35 


1.847 


1.710 


7.4% 


680 


0.061 


100 


3.005 


4.266 


41.9% 


98 x 10 3 


0.039 


1880 


6.635 


10.954 


65.1% 


34 x 10 6 



E(n 2 ,a,N) L 
3.833 
27.77 
177.5 



From table 2 we conclude that the asymptotic perturbative approximation is good 
nearly up to ni, even though the error ^(ni, a) increases significantly with decreasing a. 
Beyond ni, the two curves are numerically very different, showing that the higher-twist 
implicit in the solid curves of fig. 2, dominates. Notice also from fig. 2 that, for n fixed 
and below ni, the error decreases (slowly) with decreasing a, something expected from the 
nature of the asymptotic series. For example, 5^(35, 0.061) = 2.6% and 5^(35, 0.039) = 
2.5%. Finally, we observe from table 2 that, for all three values of a chosen, ab 2 lnni ~ 0.6, 
while ot\>2 In ?7-2 — 1 — 1.4. 



5. Summary and Discussion 

In this paper we have proposed a definition that removes the perturbative ambiguities 
of the QCD resummation formula, through a principal value prescription for integrals over 
the running coupling in the exponent E, defined in eq. (1.4). We have developed this 
approach in the context of the Drell-Yan cross section, but it should have a more general 
application. 

This resummation procedure provides an unambiguous (although arbitrary' 31 ) defini- 
tion for the resummed perturbative series for large threshold corrections in hadron-hadron 
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scattering, to which higher-twist nonperturbative corrections may, in principle, be added 
systematically. We have noted the relation of our principal value prescription to the method 
of Borel transformations. 

The exponent E(n, a) is analytically calculable, in terms of one and two-loop effects 
that contain all the large perturbative corrections in the sense discussed in ref. 2. For the 
one-loop terms of the exponent, these results are given as an infinite series of Exponential 
Integrals. This series suggests, in turn, a definition of the resummed perturbative series, 
exploiting the asymptotic approximation of these functions. The number of terms in the 
asymptotic exponent is now precisely determinable numerically, as well. 

Numerical results show that the exact leading exponent can enhance the cross section 
in a manner that is nicely approximated by its asymptotic perturbative expansion up to 
values of n such that (a s (Q 2 ) / n)b2 Inn ~ 0.6. Beyond this range higher-twist effects turn 
the exact exponent around to negative values. 

We are hopeful that this approach will be useful in organizing QCD perturbative cor- 
rections. In particular, the results for the exponent shown in fig. 2 suggest that resummed 
higher-order radiative corrections may be substantially smaller than the simple exponen- 
tiation of leading-order results. This may help to explain the relative success of one-loop 
approximations at moderate energies. In addition, we have seen that nonperturbative ef- 
fects, whose presence is implied by ambiguities in resummed perturbation theory, appear 
at order A/Q. For the normalization of the Drell-Yan cross section then, we may be in a 
situation where perturbative corrections are moderate, with relatively large, but kinemat- 
ically simple, nonperturbative corrections. Such a combination could teach us a lot about 
the interrelation of these two components of the full theory. 

Of course, the analytical result for the leading exponent, eq. (2.36), contains only an 
(arbitrary) part of the physical higher twist in dilepton production, and we have not made 
a systematic study of higher twist in this paper. A definition of the perturbative series in 
QCD is, however, necessary before higher twist can be unambiguously defined . Because 
eq. (2.36) is such a well-defined extension of the perturbative series for large corrections, 
it also contains some higher twist. This is evident numerically in the difference between 
the functions and their perturbative approximations at the edge of phase space in fig. 2. 
A complete analysis of higher twist in the normalization of the Drell-Yan cross section 



22 



is now possible, but must remain the subject of future work. One technical observation 
along these lines may be worth making here, however. In addition to the corrections at 
order (A/Q) p ,p > 1 that we have found above, we also expect corrections of the form 
1/[(1 — z)Q 2 ], associated, for instance, with the production of resonances in the final 
state. Such corrections may be equally, or even more important than those identified here, 
especially for very large Q. 

In future work we hope to explore further the numerical and phenomenological con- 
sequences of this approach, including numerical results for the next-to- leading exponent, 
and for the cross section itself. This treatment of large perturbative corrections should 
also be applicable to other inclusive hadron-hadron cross sections, including heavy-quark 

[121 

production and jet production, where gluon-gluon hard-scattering functions play an 
important role. 
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APPENDIX A 
Numerical Evaluation of the Exponent 



The first issue we will address is the number of asymptotic terms in the formula giving 
E(n, a s , N)l. As was mentioned in section 3, one definition could involve a comparison of 
each function £{mt) with its asymptotic expansion, eq. (3.2). Given that this function is 
tabulated and close to 1 at this range, this would seem a very straightforward approach. 
However, the value of N thus determined would depend on rat and one would have then to 
sum that complicated dependence over m. We took the simpler approach of using eq. (3.3) 
to define the three numbers {N[2I(t/2)], N[I(t)], N[l2(t)]} for each of the individual sums. 
These numbers of course now depend on the moment variable n, but not strongly. As we 
show in table 1, this dependence becomes more appreciable when the coupling constant 
decreases, but this has no numerical consequences for the leading asymptotic exponent in 
the perturbative regime, precisely because of the smallness of the corresponding coupling. 
The above numbers should be determined in a range of values of n that is well within the 
perturbative restriction a s lnn <C 1. The result is shown in table 1. 

Another issue is the numerical evaluation of the exact leading exponent. If we use 
the simple and straightforward analytic expression, eq. (2.36), we observe that the series 
involved* , although convergent, are alternating. Using double precision, one can not 
reliably calculate the exponent beyond a value n ~ 40, because the alternating partial 
sums exceed the 16-digit accuracy. Also, the special functions £{mt) should be calculated 
to a similar accuracy, and the corresponding numerical integrations need a lot of computer 
time. This method, for relatively low values of n, may be used as a numerical check for the 
alternative method we mentioned in section 4, namely the use of contour P, fig. 1. The 
resulting expressions for the integrals are rather complicated, from an analytical point of 
view, but run smoothly on the computer for much larger values of n. On P, the integral 
I{t), eq. (2.29), can be written as a sum of three integrals, each one resulting from twice 
the real part of the integration along one side of the contour. The result is 

I(t) = I 1 (t,r) + I 2 (t,r) + I 3 (t,r) (A.l) 
* more accurately, the finite sums involved, if we use integer n (which we will do in this case), 
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with r = n — 1, where: 



h(t, r)=t J ( ^(A 1 (r;x)B 1 (t,r; x) + Ci(r; x)D 1 (t, r; x)j 



l2M = t J^TW 



(A 2 (r ; x)B 2 (t, r; x) + C 2 (r; x)D 2 (t, r; x) ) 



l/r 

h{t,r) = -t J j^(^A 3 (r;x)B 3 (t;x) + C 3 (r;x)D 3 (t;x)^ 



and where the various integrands are given by: 



Ai(r;x)= (l+a; 2 /r 2 ] cos(r arctan(x/r)) — 1 



2 2 

B 1 (t,r;x) = (l/2)\n(l + (l/t)\n(x/r) + %/2t ) 



Ci(r;x)= (l + x 2 /r 2 j sin(r arctan(x/r)) 

/ tt/2 \ 

L>i(t,r;a;) = arctan — ; — — 

v ' ' \t + \n{x/r)) 



A2(r;x)=x(^ (1 — x) 2 + l/r 2 cos ^r arctan ~ 



l/r 



- (l/r) (l-x) 2 + l/r 2 



r/2 



sin^r arctan — j j 



S 2 (t, r; x) = (1/2) ln( [l + (l/2t) ln(x 2 + l/r 2 ) 



2 ^ _ arctan(l/xr) 
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C 2 (r;x) = (l/r)( (1 - xf + 1/r 



72 



cos r arctan 



(111 



\1 — x< 



+ X 



(1 — x) 2 + l/r 2 sin^r arctan^ 



1/r 



x 



D2{t,r;x) = arctan ( 



arctan(l/ra;) 



W + (l/2)ln(£ 2 + l/r 2 ). 



As(r; x) = x(x r cos(r7r/2) — 1) + x r sin(r7r/2) 



B 3 (t;x) = (1/2) In 1 + (1/2*) ln(l + x 



+ 



arctan x 
t 



Cs(r; x) = —(x r cos(r7r/2) — 1 — x r+1 sin(r7r/2)) 



Ds(t; x) = arctan ( 



arctan x 



Vt + (l/2)ln(l + o; 2 ). 



The integral I 2 is handled in a similar fashion. 
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FIGURE CAPTIONS 



Figure 1: The contour P used for the numerical evaluation of the leading exponent. A 
mirror-symmetric contour in the lower half-plane is understood. 

Figure 2: (a) The exact principal value leading exponent (solid curve) versus its asymp- 
totic approximation (dotted curve) as a function of n, for a(Q = 5GeV) ~ 0.075. (b) 
Similarly for a(Q = lOGeV) ~ 0.061. (c) Similarly for a(Q = 90GeV) ~ 0.039. 
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